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Abstract 



Searches for supersymmetric partners of top and bottom quarks are presented 
using data taken by the DELPHI experiment at LEP in 1997 and 1998. No devi- 
ations from standard model expectations are observed in these data sets, which 
are taken at centre-of-mass energies of 183 GeV and 189 GeV and correspond 
to integrated luminosities of 54 pb~^ and 158 pb~^. These results are used in 
combination with those obtained by DELPHI at lower centre-of-mass energies 
to exclude regions in the squark-neutralino mass plane at 95% confidence level. 
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1 Introduction 

This paper reports on a search for scalar partners of quarks (squarks) in data taken by 
DELPHI in 1997 and 1998 at centre-of-mass energies (y/s) of 183 GeV and 189 GeV. Mass 
hmits for these particles have already been published based on data taken at LEP2 0. 

Scalar partners of right- and left-handed fermions are predicted by supersymmetric 
models and, in particular, by the minimal supersymmetric extension of the standard 
model (MSSM) They could be produced pairwise via e~^e~ annihilation into Z°/7. 
Large Yukawa coupling running for the diagonal elements and important off-diagonal 
terms make the partners of heavy fermions as the most probable candidates for the 
charged lightest supersymmetric particle. As a consequence their lighter states are can- 
didates for the lightest charged supersymmetric particle. 

Throughout this paper conservation of R-parity is assumed, which implies that the 
lightest supersymmetric particle (LSP) is stable. The LSP is assumed to be the lightest 
neutralino which interacts only weakly with matter, such that events will be characterised 
by missing momentum and energy. 

In a large fraction of the MSSM parameter space sfermions are predicted to decay 
dominantly into the corresponding fermion and the lightest neutralino. Consequently in 
the search for sbottom particles only the decay into + Xi "^^is considered. For the stop 
squark, the equivalent decay into t + Xi is kinematically not allowed at LEP, and the 
decay of a stop into a bottom quark and a chargino is disfavoured in view of existing 
limits on the chargino mass 0]. The dominant two-body decay channel is thus the one 
into a charm quark and a neutralino. 

2 Detector description 

The DELPHI detector and its performance have been described in detail else- 
where I^J^; only those components relevant for the present analyses are discussed here. 
Charged particle tracks are reconstructed in the 1.2 T solenoidal magnetic field by a 
system of cylindrical tracking chambers. These are the Vertex Detector (VD), the Inner 
Detector (ID), the Time Projection Chamber (TPC) and the Outer Detector (OD). In 
addition, two planes of drift chambers aligned perpendicular to the beam axis (Forward 
Chambers A and B) track particles in the forward and backward directions, covering 
polar angles 11° < ^ < 33° and 147° < 9 < 169° with respect to the beam {z) direction. 

The VD consists of three cylindrical layers of silicon detectors, at radii 6.3 cm, 9.0 cm 
and 11.0 cm. All three layers measure coordinates in the plane transverse to the beam. 
The inner (6.3 cm) and the outer (11.0 cm) layers contain double-sided detectors to also 
measure z coordinates. The VD covers polar angles from 24° to 156°. The ID consists 
of a cylindrical drift chamber (inner radius 12 cm and outer radius 22 cm) covering 
polar angles between 15° and 165°. The TPC, the principal tracking device of DELPHI, 
consists of a cylinder of 30 cm inner radius, 122 cm outer radius and has a length of 2.7 m. 
Each end-plate has been divided into 6 sectors, with 192 sense wires used for the dE/dx 
measurement and 16 circular pad rows used for 3 dimensional space-point reconstruction. 
The OD consists of 5 layers of drift cells at radii between 192 cm and 208 cm, covering 
polar angles between 43° and 137°. 

The average momentum resolution for the charged particles in hadronic final states is 
in the range Ap/p^ ~ 0.001 to 0.01(GeV/c)~^, depending on which detectors are used in 
the track fit g]. 
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The electromagnetic calorimeters consist of the High density Projection Cham- 
ber (HPC) covering the barrel region of 40° < 9 < 140°, the Forward ElectroMagnetic 
Calorimeter (FEMC) covering 11° < ^ < 36° and 144° < 6 < 169°, and the STIC, a 
Scintillator Tile Calorimeter which extends the coverage down to 1.66° from the beam 
axis in both directions. The 40° taggers are made of single layer scintillator-lead counters 
used to veto electromagnetic particles that may be not measured in the region between 
the HPC and FEMC. The efficiency to register a photon with energy above 5 GeV at 
polar angles between 20° and 160°, measured with the LEPl data, is greater than 99% 
The hadron calorimeter (HCAL) covers 98% of the solid angle. Muons with momenta 
above 2 GeV/c penetrate the HCAL and are recorded in a set of Muon Drift Chambers. 

Decays of b-quarks are tagged using a probabilistic method based on the impact pa- 
rameters of tracks with respect to the main vertex. stands for the corresponding 
probability estimator for tracks with positive impact parameters, the sign of the impact 
parameter being defined by the jet direction. The combined probability estimator Vcomb 
includes in addition contributions from properties of reconstructed secondary vertices 0. 



3 Data samples and event generators 

Data were taken during the 1997 and 1998 LEP runs at mean centre-of-mass energies of 
183 GeV and 189 GeV, corresponding to integrated luminosities of 54 pb~^ and 158 pb~^. 

Simulated events were generated with several programs in order to evaluate the signal 
efficiency and the background contamination. All the models used JETSET 7.4 |^ for 
quark fragmentation with parameters tuned to represent DELPHI data |Q. 

Stop events were generated according to the expected differential cross-sections, us- 
ing the BASES and SPRING program packages [^. Special care was taken in the mod- 
elling of the stop hadronisation [^. Sbottom events were generated with the SUSYGEN 



program [T^. The background processes e"'"e~ — ^ qq('^7) and processes leading to four- 
fermion final states, (Z°/7)*(Z°/7)*, W"^*W~*, Wei^e, and Z°e+e~ were generated using 
PYTHIA 1^. At the generator level, the cut on the invariant mass of the virtual {Z^ /'-/)* 
in the (Z'^/7)*(Z°/7)* process was set at 2 GeV/c^, in order to be able to estimate the 
background from low mass ff pairs. The calculation of the four-fermion background was 



cross-checked using the program EXCALIBUR |[T3[, which consistently takes into account 
all amplitudes leading to a given four-fermion final state. The version of EXCALIBUR used 
does not, however, include the transverse momentum of initial state radiation. Two- 



photon interactions leading to hadronic final states were simulated using TWOGAM [|T^ and 



BDKRC [T^ for the Quark Parton Model contribution. Leptonic final states with muons 



and taus were also modelled with BDKRC. BDK |15| was used for final states with electrons 
only. 

Generated signal and background events were passed through detailed detector re- 
sponse simulation |^ and processed with the same reconstruction and analysis programs 
as the real data. The number of background events simulated is mostly several times 
larger than the number expected in the real data. 



4 Event selection 

In this section the selection to search for stop and sbottom in the decay modes ex? 
and bxi , respectively, is presented. In both cases the experimental signatures consist 
of events with two jets and missing momentum. Since event parameters, such as visible 
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energy, greatly depend on the mass difference AM between the squarlc and the LSP, 
optimised selection procedures are used for the degenerate (AM < 10 GeV/c^), and the 
non-degenerate (AM > 10 GeV/c^) mass case. The main differences between stop and 
sbottom events arise from the hadronisation, which occurs either before (t) or after (b) 
the decay of the scalar quark (in a large fraction of the MSSM parameter space the 
width of the sbottom decay into 6 + Xi is greater than the typical QCD scale so that the 
sbottom does not hadronize before it decays). These differences are visible in particular 
in the degenerate mass case. Consequently different selections are used for the stop and 
sbottom analyses in the degenerate mass case whereas the selections are identical in the 
non-degenerate mass case. 

In a first step particles are selected and clustered into jets using the Durham al- 
gorithm W^ with i/cut = 0.08. Reconstructed charged particles are required to have 



momenta above 100 MeV/c and impact parameters to the measured interaction point 
below 4 cm in the transverse plane and below 10 cm in the beam direction. Clusters in 
the calorimeters are interpreted as neutral particles if they are not associated to charged 
particles, and if their energy exceeds 100 MeV. 

In the second step of the analysis, hadronic events are selected. Only two-jet events are 
accepted. The following requirements are optimised separately for the two AM regions: 

Non-degenerate mass case: For both the stop and sbottom analyses hadronic events 
are selected by requiring at least eight charged particles, a total transverse energy^ 
greater than 15 GeV and a transverse energy of the most energetic jet greater 
than 10 GeV. These three cuts are aimed at reducing the background coming 
from two-photon processes. Forward Bhabha scattering is suppressed by requiring 
that the total energy in the FEMC is lower than 25 GeV. Z''(7) processes with a 
detected photon are reduced by requiring that the total energy in the HPC is lower 
than 40 GeV. Finally, at -^i = 183 GeV, the requirement for substantial missing 
energy is fulfilled by demanding that the quantity y/s' is lower than 170 GeV. The 
quantity is the effective centre-of-mass energy after photons radiation from 
the incoming e~^e~ beams. At ^/s =189 GeV, this requirement is replaced by the 
requirement that the polar angles of the two jets are between 20° and 160°. 

Degenerate mass case: To select hadronic events in the stop analysis the number of 
charged particles is required to be greater than five, the total charged energy has 
to be lower than 0.3^/s (in order to select events with missing energy) and the 
polar angle of the total missing momentum has to be between 15° and 165°, in 
order to reduce the background from radiative return events. The total energy in 
the FEMC and HPC has to be lower than 10 GeV and 40 GeV, respectively. The 
reduction of two-photon processes is ensured by requiring that the total transverse 



energy is greater than 5 GeV and that the quantity pu = yPtti + Ptt2 is greater 
than 5 GeV/c, where ptu is the transverse momentum of jet i with respect to 
the thrust axis projected onto the plane transverse to the beam axis. Finally, 
the most energetic charged particle is required to have a polar angle between 
30° and 150° and a momentum greater than 2 GeV/c. Similarly the polar angle 
of the most energetic neutral particle is required to be between 20° and 160°. 
At ^/s = 183 GeV, the sbottom selection at this step is similar to the stop analysis 
described above except for the requirement on pu which is replaced by requiring 
the ratio pu/Etot to be greater than 50% where Etot is the total energy of the event. 



^The transverse energy Et of a particle is defined as Et = i?| + E"^ where Ex and Ey respectively are iJcosijisinS 
and Esm4>smd . The angles </) and Q are respectively the azimuthal and polar angle of the particle. 
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At i/i = 189 GeV, the sbottom selection at this step is simphfied by removing the 
above requirement on ptt/Etot- 

After this second step and for both the non-degenerate and the degenerate mass cases, 
agreement between data and expectations from the Monte Carlo simulation describing 
standard model processes is found to be good as can be seen from Figures |l]a-c showing 
the visible mass, the charged multiplicity and the fraction of the energy for polar angles 
between 30° and 150° at y/s = 189 GeV. Figures ^-c show the total energy, the transverse 
energy and the charged multiplicity of the leading jet, for the degenerate mass case of the 
stop analysis at y/s = 189 GeV. Figures ||a-c show the visible mass, the missing transverse 
energy and the total multiplicity for the degenerate mass case of the sbottom analysis at 
= 189 GeV. 

In a third step discriminating linear functions are used in order to achieve optimum 



rejection power. They have been determined in the following way: 

Non-degenerate mass case: In this case, the same functions have been used both for 
the stop and the sbottom analysis. A first discriminating linear function has been 
determined using training samples of signal and Z'^(7) background processes. For 
the training of a second discriminating linear function, signal and WW background 
event samples have been used. In the non-degenerate mass case, these two sources 
of background processes are found to be dominant after the first and second step 
of the event selection. 

Degenerate mass case: Here the main source of background remaining after the first 
and second step of the event selection is found to be 77 events. Different functions 
have been determined for the stop and sbottom analyses using training samples 
of signal and two-photon events. 

Figure |I]d shows the discriminating function against the Z°7 background for the non- 
degenerate mass case at a/s = 189 GeV, Figure 0d and |^ show the discriminating func- 
tions for the degenerate mass domains of the stop and sbottom analyses at -^i = 189 GeV. 
For these degenerate and non-degenerate mass cases, fair agreements between data and 
expectations from Monte Carlo describing standard model processes are found. The data 
and Monte Carlo small disagreement of the discriminating function for the degenerate 
mass cases shown in Figure ^d is restricted to the negative values of this function which 
correspond to the region of the bulk of the expectations from Monte Carlo describing 
standard model processes in particular two-photon interactions leading to hadronic final 
states which are known to be difficult to modelize. This region does not correspond to the 
squark signal region. As shown by the hatched areas of Figure ^d, the positive values of 
this discriminating function correspond to the squark signal region and in this region the 
agreement between data and expectations from Monte Carlo describing standard model 
processes is very good. 

The final background reduction is performed by sequential cuts. In the non-degenerate 
mass case, one set of cuts is used to select both stop and sbottom events. It is shown, 
together with the number of events retained in data and background simulation, in Ta- 
bles and g. 

In the degenerate mass case two different selections are used for stop and sbottom. 
These are shown in Table ^ for the stop analysis at ^/s = 183 GeV and ^/s = 189 GeV 
and in Tables ^ and ^ for the sbottom analysis at y/s = 183 GeV and ^/s = 189 GeV 
respectively. 
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5 Results 

The number of candidates found and the expected background levels are shown in 
Table ^ and Table 0. There are candidates in common in the stop and sbottom analyses. 
The total background is given assuming a //or// between the degenerate mass case and 
the non-degenerate mass case for the stop and sbottom analyses. One candidate event 
from the non-degenerate mass case analysis is shown in Figure The efficiencies of the 
stop and sbottom signal selection are summarised in Figure They have been evaluated 
using 35 simulated samples at different points in the (Mq,M^o) plane, for squark masses 
between 50 and 90 GeV/c^ and neutralino masses between and 85 GeV/c^. 

No evidence for stop or sbottom production has been found in the two-body decay 
channels. Figure ^ and Figure |^ show the (Mq,M^o) regions excluded at 95% confidence 

level by the search for t — > ex? and b — > bx? decays, with the 100% branching ratio 
assumption, both for purely left-handed states (with maximum cross-section) and the 
states with minimum cross-section. We have also used the results (efficiencies, number of 
candidates and expected background) of the analyses of the data at 130 - 172 GeV [|l| in 
order to derive these exclusion regions. 

In order to estimate systematic errors coming from detector effects and modelling, the 
differences of the mean values of the observables used for the above analyses (sequential 
cuts steps and discriminating linear analyses steps) between real data and simulation are 
calculated at the level of the first step of the selection described in section ^. The difference 

6 for the mean value of the observable X between real data and simulation is used in order 
to shift X according to X + 6 and X — 6. The analyses described in section ^ then use 
the shifted observables and the differences in efficiencies and expected background with 
respect to efficiencies and expected background obtained with the unshifted observables 
are taken as systematic errors. The relative systematic errors for efficiencies are 10% 
in the non-degenerate mass case and 15% in the degenerate mass case. The systematic 
errors for the expected background are given in Table ^ and Table ^. 

Systematic errors on efficiencies coming from the modelling of the hadronization of the 
stop are estimated by switching off the hadronization of the stop. The relative systematic 
errors for efficiencies are 2% in the non degenerate mass case and 8% in the degenerate 
mass case. 

6 Conclusions 

In data samples of 54 pb~^ and 158 pb"'^ collected by the DELPHI detector at centre- 
of-mass energies of 183 GeV and 189 GeV searches are performed for events with acopla- 
nar jet pairs. The results are combined with those already obtained at centre-of-mass 
energies between 130-172 GeV . 

At 183 GeV, the search for stop and sbottom quarks, decaying into ex? and bx?, 
respectively, gives in total 3 candidates (some candidates are in common in the stop and 
sbottom analyses) well compatible with the expected background of 3.4 ± 0.5. 

At 189 GeV, the search for stop and sbottom quarks, decaying into cx? and bx?, 
respectively, gives in total 9 candidates (some are candidates are also in common in the 
stop and sbottom analyses) well compatible with the expected background of 11.6 ± 1.4. 

For the stop, a mass limit of 79 GeV/c^ is obtained for the state with minimal cross- 
section, if the mass difference between the squark and the LSP is above 15 GeV/c^ . A 
mass limit of 62 GeV/c^ is obtained for the sbottom quark under the same condition. 
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In the case of maximum cross-section, these numbers are 84 GeV/c^ for the stop and 
87 GeV/c^ for the sbottom. 
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Table 1: Fourth step of the event selection for two-body decays of stop and sbottom in the non-degenerate 
mass case at ^/s — 183 GeV. Dataigs and MCisa indicate data and Monte Carlo at ^/s — 183 GeV. DLAl 
and DLA2 denote the first and second discriminating linear analysis as explained in the text. -Pr^i^^ stands for 
the total missing momentum, Ejeti {Ejet2) denotes the energy of the (next to) leading jet, Eemjeti {Eem.jet2) 
denotes the total electromagnetic energy of the (next to) leading jet, Ojsts are the polar angles of the jets, Piao 
is the momentum of the most isolated charged particle, Othrust denotes the polar angle of the thrust axis and 
^T^h d total tranverse energy of the next to leading jet taking into account charged particles only. AM 

represents the mass difference between the squark and the LSP. The case AM > 10 GeV/c^ corresponds to the 
non-degenerate mass case. The errors on the Monte Carlo are statistical only. 
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Table 2: Fourth step of the event selection for two-body decays of stop and sbottom in the non-degenerate 
mass case at ^/s — 189 GeV. Dataisg and MCiss indicate data and Monte Carlo at -^s = 189 GeV. DLAl and 
DLA2 denote the first and second discriminating linear analysis as explained in the text. R30 {R20) denotes the 
fraction of the total energy out of the cones of 30° and 150° (20° and 160°) centered on the beam axis, p'^^^'^^a 
denotes the momentum of the leading particle; Eem2 denotes the total electromagnetic energy of the next to 
leading jet. Vcomb is a b-tagging probability as explained in the text. AM represents the mass difference between 
the squark and the LSP. The case AM > 10 GeV/c^ corresponds to the non-degenerate mass case. The errors 
on the Monte Carlo are statistical only. 
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Table 3: Fourth step of the event selection for two-body decays of stop squarks in the degenerate mass case 
at ^/s — 183 GeV and ^/s — 189 GeV. Dataisa and MCisa (Dataigg and MC189) indicate data and Monte Carlo 
at ^/s = 183 GeV (\/s = 189 GeV). R30 {R20) denotes the fraction of the total energy out of the cones of 30° 
and 150° (20° and 160°) centered on the beam axis and acoplanarity tjtrust the acoplanarity angle with respect to 
the thrust axis. For the other variables see the text as well as in Tables |l| and ^ DLA stands for discriminating 
linear analysis. AM represents the mass difference between the squark and the LSP. The case AM < lOGeV/c'^ 
corresponds to the degenerate mass case. The errors on the Monte Carlo are statistical only. 
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Table 4: Fourth step of the event selection for two-body sbottom decays in the degenerate mass case at ^/s 
= 183 GeV. Datai83 and MCisa indicate data and Monte Carlo at = 183 GeV. The variables are explained 
in the caption of Tables |l|, |^ and ^. DLA stands for discriminating linear analysis. AM represents the mass 
difference between the squark and the LSP. The case AM < lOGeV/c'^ corresponds to the degenerate mass case. 
The errors on the Monte Carlo are statistical only. 
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Table 5: Fourth step of the event selection for two-body sbottom decays in the degenerate mass case at ^/s 
= 189 GeV. Datai89 and MCisg indicate data and Monte Carlo at ^/s = 189 GeV. The variables are explained 
in the caption of Tables |l|, ^ and ^. DLA stands for discriminating linear analysis. AM represents the mass 
difference between the squark and the LSP. The case AM < lOGeV/c^ corresponds to the degenerate mass case. 
The errors on the Monte Carlo are statistical only. 
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Table 6: Number of candidates and expected background in the search for two-body decays of stop and sbottom 
when performing the //or// of the analyses in the de generate and non-degenerate mass case at ^/s = 183 GeV. 
Datai83 and MCisa indicate data and Monte Carlo at y's — 183 GeV. There are candidates in common in the 
stop and sbottom analyses. 
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Table 7: Number of candidates and expected background in the search for two-body decays of stop and sbottom 
when performing the //or// of the analyses in the degenerate and non-degenerate mass case at = 189 GeV. 
Datai89 and MCisg indicate data and Monte Carlo at y/s — 189 GeV. There are candidates in common in the 
stop and sbottom analyses. 
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Figure 1: a) the visible mass, b) the charged multiplicity, c) the fraction of the energy in the polar angle 
interval [30°, 150°] and d) the discriminating function against the background (as described in the text) for 
the non-degenerate mass case concerning both stop and sbottom analysis. The dots with error bars show the 
data while the clear histogram is the SM prediction. Each hatched area shows the stop signal for stop masses of 
70 GeV/c^ , 80 GeV/c^ and 90 GeV/c^ with AM > 10 GeV/c^ (with a normalization factor to the luminosity 
in the range 8 to 90) where AAf represents the mass difference between the squark and the LSP. The case 
AAf > 10 GeY /(? corresponds to the non-degenerate mass case. 
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Figure 2: a) the total energy, b) the transverse energy, c) the charged muhipUcity of the leading jet and d) 
the discriminating function (as described in the text) for the degenerate mass case of the stop analysis. The dots 
with error bars show the data while the clear histogram is the SM prediction. Each hatched area shows the stop 
signal for stop masses of 70 GeV/c^ , 80 GeV/c^ and 90 GeV/c^ with AM < 10 GeV/c^ (with a normahzation 
factor to the luminosity in the range 8 to 90) where AA'I represents the mass difference between the squark and 
the LSP. The case AM < 10 GeV/c^ corresponds to the degenerate mass case. 



14 



V s =189 GeV AM < 10 Ge\/c 




25 50 75 100 2 4 6 8 10 

2 

Visible mass (GeV/c ) Missing transverse energy (GeV) 

F 900 r 




10 20 30 40 50 -4 -2 2 4 

Total multiplicity Discriminating function 



Figure 3: a) the visible mass, b) the missing transverse energy, c) the total multiplicity and d) the discriminating 
function (as described in the text) for the degenerate mass case of the sbottom analysis. Each hatched area shows 
the sbottom signal for sbottom masses of 50 GeV/c^ , 60 GeV/c^ , 70 GeV/c^ , 80 GeV/c^ and 90 GeV/c^ with 
AM < 10 GeV/c^ (with a normalization factor to the luminosity in the range 5 to 100) where AM represents 
the mass difference between the squark and the LSP. The case AA{ < 10 GeV /c^ corresponds to the degenerate 
mass case. 
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Figure 4: View of one candidate event from the non-degenerate mass case in the transverse plane. The 
corresponding total energy is 57.3 GeV, the charged multiplicity is found to be 27, the total visible mass is 43.3 
GeV/c^, the polar angle of the missing momentum is 74.8 degrees and the polar angle of the two jets are 86.5 
degrees and 125 degrees respectively. 



16 




10 20 30 40 50 60 70 80 90 

LSP mass (Ge\/c) 



Figure 5: Efficiencies for tlie a) sbottom and b) stop selection in the search for two-body decays as function 
of tire LSP mass for various sbottom and stop masses. The sbottom and stop masses are indicated on the plots 
in units of GeV/c^. 
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Figure 6: Exclusion domains at 95% confidence level in the (t,Xi) mass plane assuming 100% branching ratio 
into cx° for pure left-handed state = rad) and for the minimum cross-section {9 — 0.98 rad) corresponding 
to the decoupling of the stop from the Z boson. The limits are obtained combining data at ^/s — 130 - 189 GeV. 
The shaded areas have been excluded by LEPl @ and CDF jl|]. 
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Fi gure 7; Exclusion domains at 95% confidence level in the (fe,Xi) niass plane assuming 100% branching ratio 
into bxi for pure left-handed state (6 = rad) and for the minimum cross-section (6 = 1.17 rad) corresponding 
to the decoupling of the sbottom from the Z boson. The limits are obtained combining data at ^/s = 130 - 189 
GeV. 



